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We study the influence of magnetic scattering on the Josephson critical current, R, in a 
superconductor/ferromagnetic metal/superconductor (SFS) junction by a tunneling Hamil- 
tonian approach. An analytical formula of R is given in the fourth order perturbation theory 
as regards the tunneling matrix element. The R exhibits the damped oscillatory dependence 
on the thickness of the ferromagnetic metal, d, and shows the transition between 0- and 
TT- states with d. When the superconducting transition temperature is comparable to the fer- 
romagnetic Curie temperature, the period of oscillation is obviously changed by increasing 
temperature, T, due to the magnetic scattering, which induces the O-tt transition with T. 
The magnetic scattering provides rich variety of Josephson effect in the SFS junction. Our 
results present an appropriate condition of a superconductor and a ferromagnetic metal to 
control the 0- and the tt- states. 



1. Introduction 

The Josephson effect is characterized by the current through a thin insulating layer without 
voltage-drop.^ It is a macroscopic quantum- mechanical phenomenon to preserve the phase 
coherence between two superconductors. A relation between the current, I, and the phase di 
fference, (p, is given by / = Ic sin (p. A finite voltage-drop appears for larger values of I than the 
Josephson critical current, R, which decreases monotonically with temperature, T.^ Similar 
phenomenon is observed in a weak link through a metallic layer due to the proximity effect.^ In 
a superconductor /normal metal/superconductor (SNS) junction, R decreases monotonically 
with the thickness of metallic layer as well as with T.^ 

There are growing interest in a superconductor/ferromagnetic metal/superconductor 
(SFS) junction, in which R shows a cusp as a function of thickness of ferromagnetic layer, 
d, and/or Such a non-monotonous behavior is in marked contrast to R in the SNS junc- 

tion. Since the cusp originates from taking an absolute value of the oscillation with alternating 
sign, the current-phase relation (CPR) is considered to be shifted by tt in a certain range of d 
and/or T in the junction from that in the SNS junction. It is called Tr-state, while the conven- 
tional Josephson junction has 0-state. The rr-state is experimentally observed by measuring 
the CPR or a spontaneous magnetic flux in a superconducting ring including the SFS junc- 
tion. The CPR is obtained by minimizing the free energy, F, as, I = 2e{dF/dp>) = R sin ip. 
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Then, the tt - state {0-state) has a negative (positive) coupling constant as /c < 0 (/c > 0). 

Such a negative coupling constant was originally found in the Josephson effect with a spin- flip 

21 - 2 “^ 

process. 

The -K-state in the SFS junction, that was first predicted by Buzdin et al?^ is a novel 
phenomenon obtained by combining a superconductor (SC) with a ferromagnetic metal (FM), 
although these orders generally compete with each other. Recently, many theoretical studies 
are made on the tt - state in the SFS junction. So far, semiclassical approaches like 
the Eilenberger equation^^ and the Usadel equation^^ are adopted by many authors. 

The Bogoliubov-de Gennes equation provides another viewpoint on the n-state in terms of 
Andreev bound state. A key mechanism is quite similar to that of Fulde-Ferrell-Larkin- 
Ovchinnikov (FFLO) state. Cooper pairs penetrating in the FM acquire a finite center- 
of-mass momentum proportional to the magnetic exchange splitting, fiex(T), between the up- 
and down-spin bands. As a result, the pair correlation oscillates in the FM as a function of 
distance from the interface. This oscillation provides the vr-phase shift depending on d. For 
example, if d is a half of the period of the oscillation, the CPR is shifted by tt from that of 
a normal Josephson junction. Hence, the transition between the 0- and vr states occurs by 
changing d. However, this period is not determined only by he^{T) but also by randomness 
and spin-dependent scattering in the In particular, the latter process plays an 

important role in the T-induced O-vr transition. It is noted that the n-state in the SFS junction 
remains alive even in a dirty FM, whereas the FFLO state is fragile against to a tiny amount 
of randomness. 

In this paper, an analytical formula of Ic in the SFS junction is presented in the fourth 
order perturbation theory by use of tunneling Hamiltonian^®^^® and Green’s function with 
path integral framework. We study the influence of magnetic scattering on the d- and T- 
dependences of 1^ in which the transition between the 0- and vr-states occurs. It is found 
that the magnetic scattering is important to observe the O-vr transition, in particular, with 
T. The tunneling Hamiltonian method makes it possible to treat systematically the magnetic 
scattering from a microscopic viewpoint. 

As an application of the -K-state in the Josephson junction, several authors have proposed 
a qubit, that is a quantum coherent two-level system and is utilized in a process of quantum 
computing.^^’^®’^^ In general, however, the qubit made with the Josephson junction has a 
short decoherence time due to several couplings to internal degrees of freedom and/or to 
the environment. The SFS junction has more degrees of freedom like a magnetic scattering 
than the conventional Josephson junction to break the coherence of Cooper pair. Our results 
provide conditions to control the 0- and vr-states and would show a route to realize the solid- 
state qubit. 

The rest of this paper is organized as follows. In §2, we briefly summarize our model 
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Hamiltonian and present a formulation of the Ic in the SFS junction. In §3, the d-dependence 
of Ic is calculated for several conditions. A role of magnetic scattering is clarified by showing 
T-dependence of Ic for several cases. Summary and discussions are given in §4. Below, h = 1 
and A:b = 1 are used in the equations. In order to obtain an actual value, h = 6.6 x 10 “^® eV-s 
and k-Q = 8.6 x 10“® eV-K“^ are used. 

2. Tunneling Hamiltonian approach for Josephson critical current 

2.1 Model Hamiltonian 

We consider s-wave SC’s and a FM. The mean- field Hamiltonian, i7, is given by 



H 


= FFsc + FFfm + FFt, 


(1) 


Hsc 


= Hr + Hr, 


(2) 


Hl 


kha 






+ ^Ae‘^L +h-c., 
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tjR 
11 ^ 




(15) 



In eq. (2), (F^r) is the Hamiltonian of left (right) SC, where i^kn„) is the kinetic energy 

of electrons with wavenumber feucr (^Ro-) and /x is the chemical potential. The electron effective 
mass is denoted by m. The superconducting gap and the phase variable are denoted by A and 
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(^L (v^r)) respectively. The annihilation and creation operators of electron with wavenumber 
vector fco- and spin a are denoted by and cl^, respectively. Equation (7) describes the ki- 
netic energy of electrons in the FM with the magnetic exchange energy, he^(T) = (Jh/2)(5'^), 
where (S^) is determined by the self-consistent equation as, (S^) = (1/2) tanh[2(TFM/T)(5^)], 
and Jh is a coupling constant between an electron and a localized moment in the FM with N 
being the number of atomic sites. The Curie temperature, Tfm; is determined by the ferro- 
magnetic coupling constant, Jy < 0, in eq. (9). For the ferromagnetically ordered moments in 
eq. (9), we adopted a mean-field Hamiltonian given by Jy Y2{ij)i^i) ' ^ where {i,j) indicates 

the sum of nearest neighbor sites, and Sf denotes the z component of the localized spin at 
the i-th site. In eq. (10), a non-magnetic impurity potential, u{q), in the FM is averaged in 
the order of UiU^, where rii is the density of non-magnetic impurity and u{q) is assumed to be 
independent of q. We adopt eq.(ll) to describe a magnetic scattering, SS^ = — (Sg), in the 

FM, where q means the wavenumber. The tunneling Hamiltonian is given by in eq.(13), 
where (r*/?) is the position of the interface between the left (right) SC and the FM. 

The current flows in order to minimize the coupling energy between two SC’s, and is given 



by 



dF 

I = 2e— = /csmv9, 

dip 



(16) 



where F is the Free energy, and (p = ipi^ — is a phase difference between two SC’s. Hence, 
it is F as a function of p> that we need to calculate. Details are summarized in Appendix A. 
2.2 Clean system 

In the clean system, the fourth order term of F as regards t is shown in Fig. 1. Detailed 
calculations are shown in Appendix A. For hex{T)/p, LVn/ p ^ 1, the analytical form of Ic is 
given by, 



IrRn = A" 



rE 



1 

u;2+A2* 



^-2\u}n\d/v-p 



cos 



2/lex (T) 



UF 



(17) 



where Rq^ = {mV /2 is a constant determined by the material and the interfaces. 
The density of states at the Fermi energy is denoted by A^f> and V is the volume of FM. 
The Ic is plotted in Fig. 2 for uf=2.5x10^ m/s, /io=0.36 eV, A=1.5 meV, T=4 K, and 
Tfm/Tsc=10. ho is defined as, /iex(T=0K) = ho- The ferromagnetic and the superconducting 
transition temperatures are denoted by Tfm and Esc, respectively. When d increases from 
d=0, in which the ground state should be the 0-state, d=nvp / 4hex{T) is the thickness for 
the first minimum of R- It is found that only the ratio of hex{T) and uf provides the period 
of oscillation in the R-d curve. If Ni or Co are chosen as the FM, Tfm is much larger than 
Tsc, i-e., Tfm/Tsc 1, and /lex(T) is almost equal to ho below Tsc- Therefore, the period of 
oscillation is estimated as, 7 tuf//io= 1-4 nm, for uf=2.5x 10^ m/s, /io=0.36 eV. This result is 
consistent with experimental results using clean FM.^^^^® Although another combination of 
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uf=2.5x10® m/s and /iq=3.6 eV also provides the same period as i:vp /hQ=lA nm, this value, 
i.e., /io=3.6 eV, is considered to be too large for hex{T). Then, we will use up=2.5xl0^ m/s, 
below. 

The Ic decays with T and/or d, whose dependences are determined after the summation 
as regards ujn in eq. (17). It can be done as, 

^-2\u)n\d/v-p 
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where = uf/27tT. See also Appendix B. By considering eq. (18), the T-dependence of 
eq. (17) is shown in Fig. 3 for d=5 nm and 100 nm. If the thermal length, is much smaller 
than the thickness, d, i.e., ^ d, eq. (18) is estimated as, 
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^T/ Ct 

This result in eq. (19) is shown by the broken line in Fig. 3, which is obtained by replacing 
the sum of LOn in eq. (17) with uJn=o = vrT as well. In the high temperature region defined by 
d/^T 3> 1, the Cooper pairs thermally lose their coherence in the FM and then Ic exponentially 
decays as exp(— d/c^x)- On the other hand, the solid curve in Fig. 3, which is the Ic-T curve 
for d/^T 1) deviates from the exponential behavior and looks like the Ic-T curve in the 
SIS junction given by Ic = d?j^^Atanh(A/2T). denotes normal resistance of insulating 
barrier.^ This similarity means that the Cooper pairs can go through the FM before losing 
their coherence, if the FM layer is thinner than the thermal length. 

We can show R for d/^x ^ 1 at T=0 K as, 
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( 20 ) 

( 21 ) 



where .^o = uf/2A. It is noted that R in the clean system at T=0 K decays with d in the 
power law as, Co/d. 

2. 3 Dirty FM with magnetic scattering 

In a dirty FM including the magnetic- and non-magnetic scatterings, the fourth order 
term of F is calculated as shown in Fig. 4. The diffusive motion of electrons in the FM is 
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indicated by F, that is given by a sequential sum of the non- magnetic impurity and magnetic 
scatterings as shown in Fig. 5. The double solid lines represent the Green function in the FM 
with the self-energy, S (wn), in the Born approximation given by 

i 



S (oJn) = - 

1 



' imp 



2 \Ti 
= 2nNpniU 



1 1 
— + 



imp 
2 



r.{T) 



sgn (uJn) , 



1 



= 27tNf(^] {5S^5S^) = — {5S^5S^). 

' Ts ( T ) V 2 / TsO 



( 22 ) 

(23) 

(24) 



The relaxation times by the non-magnetic and magnetic scatterings are denoted by Timp and 
Ts(T), respectively. In this study, the spin-spin correlation function is approximated by that 
of Ising-type mean-field Hamiltonian as, 



(SS^SS^) = 1 



^ cosh" 



1 -1 



(5.) - 1 



(25) 



_Tfm V T 

By taking the dirty limit of hex{T)T, vp\Q\r <C 1, one can obtain the following equation 
as, 



F = -2ttNp 



JQd 



Q 



DQ^ + 2 \cOn\+ 2/rs(T) ± i2he.{T) ’ 
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— + 



-1 



(26) 



(27) 



‘ imp rs(T), 

where D is the diffusion constant in the FM. It should be noted that D contains Ts(T) as well 
as Timp, while D has been so far assumed to be constant by other authors. We find that the 
diffusive motion of electron can be changed by T due to the T-dependence of Ts{T). We will 
discuss this point in the next section. 

Finally, the Ic in the dirty FM is given by 

1 



IcRd = A^T 






cut; + 



exp 



d\ f d 




(28) 
(29) 

y(cun + l/rs(T))2 + /lex(T )2 + (cu„ + 1 /Ts(T))J ’ 

where = t^NpTr‘^V/2Dd is assumed to be a constant determined by the material and the 
interfaces. The summation in eq. (28) is numerically carried out, and results are shown in the 
next section. 

3. d- and T-dependences of 0- and 7r-states 

We have shown in eq. (17) that the R in the clean system oscillates as a function of d. 
The period is determined by 2he^{T)/vp, and the decay length is given by for d/^T 3> 1. 
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On the other hand, in the dirty FM with magnetic scattering, the period, and the decay 
length, depend on /lex(F), D, and Ts(T), as shown in eqs. (29) and (30). The |/c|-<^ curves 
are plotted in Fig. 6 for (a) Tfm/Tsc= 10, ho=30 meV, Timp=2xl0“^^ s, Tso=10“^^ s, (b) 
Tfm/Tsc= 1, ho=30 meV, Timp=2x 10“^"^ s, T50=10“^^ s, and (c) Tfm/Tsc= 1, ho=100 meV, 
Tmp=2xl0“^^ s, s. The other parameters are set as, A=1.5 meV, Tsc=10 K, and 

uf=2.5x 10^ m/s. Although the ir-state is characterized by a negative value of /c, |/c| is useful 
to compare our results with experimental ones. 

In all cases, |/c| shows the damped oscillatory behavior as a function of d, and the 0 -tt 
transition occurs. For ho=30 meV at T/Tsc=0.2 in Figs. 6 (a) and (b), the first minimum 
appears around d=6.3 nm, which can be considered the case of ferromagnetic metal-alloy, 
PdNi.® If one assumes a larger value of ho as 100 meV, the period of oscillation becomes 
small. However, by taking a smaller value of Timp as 2x10“^^ s, the first minimum appears 
around 7.5 nm as shown in Fig. 6 (c), since is determined by both he^(T) and D including 

For Tfm/Tsc= 10 in Fig- 6 (a), the period of oscillation does not change with T. On the 
other hand, in the case of Tfm/Tsc= 1) the period is strongly elongated or shortened with T 
as shown in Figs. 6 (b) and (c). The difference between the two cases, i.e., Tfm/Tsc= 10 and 
1, originates from the T-dependence of Ts{T) and hex{T) in eq.(29). For Tfm/Tsc=1) Ts{T)~^ 
{he^(T)) exponentially increases (decreases) with T as shown in Fig. 7, while these are almost 
constant below Tsc for Tfm/Tsc= 10. Such T-dependences are attributed to the magnetic 
scattering. Several authors have pointed out that the magnetic relaxation time elongates the 
period of Ic and changes the ratio of the decay length and the period.^®’ Our result 
includes their one, which corresponds to the case of Fig. 6 (b). Some experimental studies 
have reported that the cusp in the \Ic\-d curve is shifted toward a larger value of d with 
increasing 16 g in which Tfm/Tsc= 1 and Timp 7 ^ 50 = 0 . 1 , is consistent with 

these experimental results. On the other hand, in a certain case like Timp/'Tso=l, the period 
of oscillation decreases with T, and the minimum of |/c| is shifted toward zero as shown in 
Fig. 6 (c). Although this case is not yet observed in an experiment, it would be realized by 
controlling the randomness and the magnetism. 

What does make this difference between Figs. 6 (b) and (c)? Since D is composed of Timp 
and Ts{T) in eq. (27), the T-dependence of D changes depending on the ratio, rimp/'rs(T). For 
example, if 1 /ts(T) <C l/rimp, D can be approximated as a constant, while D exponentially 
decreases with T for 1 /ts(T) ^ l/uimp- This fact means that the ratio of non-magnetic- 
and magnetic relaxation times, Timp/'rs(T), is also an important parameter to control the T- 
dependence of |/c|- As we have shown by the perturbative calculation that Ts(T) is included 
in the denominator of and in D in eq. (29), the period of oscillation, shows various 
behaviors with T and rimp/'rs(T). In Fig. 8, is plotted for (a) ho=30 meV, rimp=2xl0“^^ 
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s, and (b) /io=100 meV, Timp=2xl0“^^ s. Other parameters are set to uf=2.5x 10^ m/s, 
2 fm/7sc =1) and Tso=10“^^ s. In Fig. 8 (a), which provides Fig. 6 (b), increases with T, 
while it starts to decrease for Timp > 0.5x10“^^ s in Fig. 8 (b). 

By setting d near the thickness of /c ~ 0, the O-vr transition occurs by changing T for 
^fm/7sc= 1) while it is not found for Tfm/2sc= 10. In Fig. 9, the T-dependences of |/c| are 
shown for (a) d=6.5 nm and (b) d=7.5 nm. In both cases, the ground state is the Tr-state. The 
solid- and broken lines show |/c| for Tfm/Tsc= 1 and for Tfm/Tsc= 10) respectively. In the 
latter case, |/d monotonously decreases with T as found in the SIS junction, since the period 
of oscillation, does not change with T. On the other hand, when Tfm is comparable to 
\Ic\ shows the cusp in the curve, since is strongly elongated or in some cases shorten 
with T as shown in Figs. 8 (b) and (c). As a result, Ic changes its sign from negative to positive 
by increasing T and the 0 -tt transition occurs. The T-induced O-vr transition is observed in 
some experiments,®A0Ti,i3, 14, 16 can be attributed to the magnetic scattering leading 

to the T-dependence of with Tfm ~ Tsc- 

We notice that |Ic| in higher T is enhanced compared to that in lower T in Fig. 9 (a), 
and vice versa in Fig. 9 (b). If we remembered the |/c|-T curve in the SIS- and/or the SNS 
junctions. Fig. 9 (a) would seem to be strange. However, the thickness in Fig. 9 (a) is closer 
to the cusp in the |Ic|-d curve than that in Fig. 9 (b). Thus, due to the large shift of cusp 
with T in the \Ic\-d curve. Fig. 9 (a) also becomes possible. Actually, both behaviors shown 
in Figs. 9 (a) and (b) are observed in the case of Cu 52 Ni 4 g with Tfm < 20 K and Nb with 
Tsc = 9.23 K by Sellier et al.^^ 

The ratios of the period and the decay length with cu„=o = vrT is given by 

'7rT + Ts(T)-i' 



^1 = 1+2 



+ 2 



heAT) 

7rT + Ts(T)-i^' 



1 + 



7rT + Ts(T) 

heAT) 



-1 



(31) 



1-14 



heAT) 

which is plotted in Fig. 10; (a) by the solid line for Tfm/Tsc=10) /iq= 30 meV, Timp=2x 10“ s, 
Tso=10“^^ s, (b) by the broken line for Tfm/Tsc=1) ho=30 meV, Timp=2x 10“^^ s, Tso=10“^^ s, 
and (c) by the dotted broken line for Tfm/Tsc=1) /io= 100 meV, Timp=2xl0“^^ s, T50=10“^^ 
s. In all cases, ^-/^+ linearly increases with T around T/Tsc ~ 0, since heAT) and Ts(T) 
do not change with T and only oJn=o = tiT contributes to the T-dependence of ^_/^+. For 
Tfm/Tsc= 10, ^-/^+ linearly increases upto Tsc due to ttT > 1/ts(T). On the other hand, in 
the case of Tfm/Tsc=1) i-/^+ are exponentially enhanced below Tsc due to T-dependence of 
heAT) and Ts(T) as shown in Fig. 7. Some authors have claimed that is approximately 

equal to 1 for heAT) 1+ and some experimental studies have reported that ^_/^+^l 

in the Josephson junction with dirty Our results for heAT) ^ T, 1/ts(T) also show 

^_/^_i_i^l, which is consistent with those theoretical and experimental results. Moreover, we 
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have shown that C-/C+ can be enhanced by T due to the magnetic scatterings, which provides 
new possibilities and rich variety for the Josephson effect in the SFS junction as shown in this 
paper. 

4. Summary and discussions 

In this paper, we have studied the Josephson critical current, Ic, in the superconduc- 
tor/ferromagnetic metal/superconductor (SFS) junction by a tunneling Hamiltonian method. 
The analytical formula of Ic is presented by the Green’s function with path integral framework 
in the fourth order perturbation theory as regards the tunneling matrix element. The influ- 
ences of magnetic scattering on the 1^ are discussed by calculating 1^ as a function of thickness 
of ferromagnetic metal (FM), d, and temperature, T. The Ic exhibits the damped oscillatory 
dependence on the thickness, and shows a transition between 0- and n-states. The oscillation 
period in the clean system is determined by the magnetic exchange splitting between the up- 
and down-spin bands in the FM, he^(T). In the dirty FM with magnetic scattering, the period 
depends on hex{T), diffusion constant, and magnetic relaxation time, Ts(T). In particular, we 
found that Ts{T) plays an important role in the temperature-induced O-vr transition. When the 
superconducting transition temperature is comparable to the Curie temperature, it is found 
that the period of this oscillation is obviously changed by increasing T due to the magnetic 
scattering. 

In our approach, the interfaces between FM and superconductors (SC’s) are tunnel-like 
and their order parameters sharply drop at the interface. On the other hand, the spatial 
variation of order parameters on the SC side near the interfaces becomes important for low 
tunnel-barrier or junction with a metallic interface. Since our method does not need any 
special boundary condition at an interface, such a variation can be included by taking higher 
order terms about the tunneling matrix element and self-consistent calculations. However, we 
will leave it in a future issue. 

Recent progress in the growth techniques for thin films provides various types of mag- 
netic Josephson junctions like a SC/anti-ferromagnet(AFM)/SC junction,®^ and a SC/half- 
metallic ferromagnet/SC junction. Moreover, some of multilayered high-Tc cuprates like 
HgBa 2 Ca 4 Cu 50 y inherently contain the magnetic Josephson junction, in which AFM sepa- 
rates two Related to the high-Tc cuprates, the Josephson junction of d-wave SC 

has been studied by some authors. Researches on these magnetic Josephson junctions will 
open up a new possibility of interplay between the Josephson current and magnetism. Our 
method is applicable to such various types of magnetic Josephson junctions.®^’ 

Finally, although the spin scattering is partly included in this study, the spin wave ex- 
citation is not considered. Such massless excitations will generate additional dissipation in 
the Josephson current, which is observed as a tiny voltage drop below Ic in the SC/normal 
metal/SC (SNS) junction.^^’®^’®^ The dissipation in SFS junctions is an interesting research 
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topic in order to realize a solid state qubit. 
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Appendix A; Perturbative calculation of 1^ in path integral framework 

A.l Basic formula 

In the path integral framework, the partition function is given by 



where 



Z = TV[e“'’"| = = J , 



S= T. 






= 



V'LT(^LTdWn) 
V'FMt(^FMt) i<^n) 



Go' = 



0 Gl 



0 ^?R^ 



Each element of Gq ^ is the Green’s function in the SC, 



^L(R)cr “ 2m ^ J 

-AL(R)e''^FH)-3^^^ 



and that in the EM, 



G'FMa = -i^^nO-Q “ M “ Crhe^{T) j (T3, (A-6) 

where aa (a=0, 1, 2, 3) is the Pauli matrix. In this framework, the tunneling Hamiltonian is 
given in the following matrix form as, 



T = 



0 Ti 2 0 

T 21 0 T 23 

0 T32 0 



Ti 2 = -fexp[i(fcFM<7 - fcL<7)’’i]o-3 = 721, 



723 = -texp[i(fcR^ - fcFMa)»’R]o-3 = T3* 
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The phase factor in the tunneling matrix element originates from the position vector of the 
interface, since a flat interface between SC and FM is written as, t6{r — tl). After integrating 
T out in eq. (A-3), the free energy is given by 

F = -Trln[-Go^ +T], (ATO) 



In the fourth order perturbation theory about t, F is approximated to be, 

F ~ Trln[Go^] 

+ ^Tr[GoTGoT] + ^Tr[GoTGoTGoTGoT]. (ATI) 

The first and second terms do not contribute to the Josephson coupling energy, whose leading 
term is the third one in eq. (ATI). 

A. 2 Clean system 

In the clean system, the forth order term is given by 

Feff =2t"^rRe[ Q-KkFMa-kFM-aHrR-rL) 






^ fRa{k^cr,i^n)9py[cr^kFMcTj '^a)n) 

^ 5fM— CT) ~^^n) fpo-ikLaA^n)]i 



where each Green’s function is given by 

/L(R)cr(^L(R)o-T<^n) = 
fl'FM(T(^FMo-T^n) = ~T 



AL(R)e'^TR) 



+ '^L(R)n- + Aj, 



L(R) 



I 



A'fMo-^ kpMai i^^n) — 



i<^n — ^FMo- ’ 

1 



iWn + ■^FMo 



^L(R) — 



^L(R)o 

2m 



IJ', 



«PM = 

Note that fepMo- and kpM-a are independent of each other. 
Finally, we obtain the following formula as, 



xRe 



E 



gSgn((.t)„ )id.y/ 2m{fi+hex {T)+iijj„ ) 

c^^ + A2 (T) 



^-sgn{ijj„)id^y 2m{fi+heFF+'^‘^n) ^iip 

^UA2 (T) 



(AT2) 



(ATS) 

(AT4) 

(ATS) 

(AT6) 

(AT7) 



(ATS) 
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For hex{T)/fi, ^ 1, eq. (A-18) is approximates as 

(mYNp)^ ^ 






-FA\T) 



rV - 



,-2\ u ,„\ d / vF \ PQg f ‘^ he ^{ T)d 



cvl + A^T) 



Vp 



cos((p), 



(A-19) 



where V = S • d and S is the cross section between SC and FM. 

A. 3 Disordered system with magnetic scattering 

In the dirty FM with the magnetic scattering, the forth order term is given by 
Feff = 2t^iRe[ ^ ^--^QirR-rD 

{k},\LOn,a 

xfRa{kRcr, i^^n)GFMa{kFMcr, ii^n) 



GfM— (j( ^FM— (7) i^n)/L(j (^L ct) i^n)F(Q, iWji)] ) 



(A-20) 



where Q = kpMa ~ kpM-a in the above notation. The Green’s function in the dirty FM is 
given by 

1 



GfMct (^( 7 ) i^n) — 



Cfcc Zj(J 

where the self-energy, So-(fc,io;n), is approximated as, 

S^(fc^,iwn) = -i(^;^ + ^^l^^ )sgn(cun), 

rs(T)-i = 27r{^f^NF{6S^6S^). 

The diffusive motion of electrons in the FM is described by F defined as, 
F(Q,iu;„) = ^A™(Q,icu,), 

m 

Q - iQ { rR - rL ) 



(A-21) 

(A-22) 

(A-23) 



-2Nf7t 



K(^Q , itOn) 



Q 



|2 , 

= “I”! + 



DQ^ + 2 |wn| + 2 /ts(T) ± i2hex(T) 



= F, 



(A-24) 






X G-i{k-^,i‘^n)Gi{—ki + Q, —iujn), 

k(j 

Amp “ A(^) 



Fimp T Ts{T) 



(1-2 \uJn\ T - DQ^t =f i2/iex(F)r) . 



(A-25) 

(A-26) 
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In eqs (A-24) and (A-26), we take the dirty limit as, oJnT, hex{T)T, vy\Q\t <C 1. Finally, we 
obtain the free energy in the dirty FM as, 

t^N^TT^V 



FeS — - 



ADd 
xRe[T^ 



■A^(T) 
1 



u;2+A2(r) 



(e-“+'^+e-““'^)e‘^]. 



= 



2(|w„| + 1 /ts(T) ± ihex(F)) 



D 



Appendix B: Summation of Matsubara frequency in 

A periodic function given by, 



(A-27) 

(A-28) 



f{t) = -/(t + (2n + l)/r), 

= /(t + 2n/T), 



can be expanded as, 

f{t) = T^e--"*, 

n 

(jJn = (2n + l)7rT. 

By taking account of eqs. (B-1) and (B-2), an alternating pulse function defined by 






is transformed as. 



= J — e“^*(5i/'r(f) — 6{uj — ujn)- 



Equation (B-5) is useful to sum up the Matsubara frequency in eq. (17) as, 

^-2\oj„\d/vp 






a;, 



2 + A2 



^-2\u}n\d/v-p 



2 + A2 



E 

n 



-5{ui - UJn) 



^—2\u)\d/vY-\-i^t 



UJ 



2 + A2 



oo ^-2u)d/vY-\-i^t ^—2u}d/v-p—iu}t 

dcu 



1 [ ^ 
A / ^ 



poo 

di/T{t) / dxsi 
JO 



smx 



X 



1 



+ 



Cj2+A2 



1 



[x + {2d/vY — it)A X + {2d/vY + h)Aj ’ 
1 r dt 



A 






dx sin X 



(B-1) 

(B-2) 

(B-3) 

(B-4) 

(B-5) 



(B-6) 

(B-7) 

(B-8) 

(B-9) 
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x + d/Co 

(x + d/Co)'^ + (tA)2 



1 

2^ 



dx sin X 



xE(-i) 



X + d/^o 

(x + d/Co)2 + n2(A/T)2 



(B-10) 



(B-11) 



1 



27tA 



/-OO 


'ttT , , 


f ttTx 


d \ 


/ dx sin X 
/o 


——cosech 
A ' 


\ A 





(B-12) 



where = vpI27tT and = up/2A. The following relation is used to obtain the forth 
equality in the above equation, 



duj 



^-piv 



u;2 + 1 



= / dx 



sinx 
X + p 



= — cos(p)si(p) + sin(p)ci(p). 

In the limit of d/^p 3> 1, eq. (B-12) is approximated as, 

1 /■“ . 

dx smx 



27tA 

1 

2^ 





^ ttTx 


d \ 


——cosech 

A ' 


K A 







'2ttT / 


^ ttTx 


V 


sinx X 


i A 


. A 


^T/_ 



oc exp { — ) for d/^p » 1. 
ST, 



On the other hand, in the limit of T ^ 0 K, 



lim cosech 
T^OK 



ttT 



X + 



2dA\ 



VF 



A 



ttTx + d/^o' 



where = xf/ 2A. Then, 



lim T 



oc 



T^O 

1 

2^ 

1 

2^ 

1 

d/^o 



^-2\uiri\d/vY 

c^^ + A2 ’ 
sinx 



dx 



X + d/io 



— cos 



SI 



+ sin 



^0 



Cl 



(B-13) 

(B-14) 

(B-15) 

(B-16) 

(B-17) 

(B-18) 

(B-19) 
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Fig. 1. Schematic figure of SFS junction and the fourth order diagram for F . The circles with up- and 
down arrows in SC indicate a Cooper pair, and the circle with up (down) arrow in FM indicate 
the up-spin electron (down-spin electron). A tunneling matrix element is denoted by t. 




Fig. 2. Ic in the clean system for uf=2.5x 10® m/s, /iq= 0.36 eV, A=1.5 meV, T=4 K, and 
7fm/7sc=10. 
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Fig. 3. The critical current Ic as a function of T is plotted for wf=2.5x 10® m/s and Tsc=10 K 
with d=5 nm by the solid line and with d=100 nm by the broken line. Here, is greater than 
vf/2ttTsc —30 nm. Hence, the solid line always satisfies d/^T <1 and becomes like the Ic~T curve 
in the SIS junction, while the broken line seems to be that in the SNS junction due to c?/^t >1 
for T/Tsc^O.3. 




Fig. 4. (a) The fourth order diagram contributing to the Josephson current in the dirty FM. t is a 

tunneling matrix element. The double solid lines represent the Green function with the self-energy, 
S (wn)) in the Born approximation, (b) The self-energy by non-magnetic and magnetic scatterings 
in the Born approximation. 
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Fig. 5. Diagram of diffusive electrons in the FM by non-magnetic and magnetic scatterings. 
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Fig. 6. The absolute values of Josephson critical current as a function of d for ■(;f=2.5x 10^ m/s, 
and (a) Tfm/Tsc=10, ho=30 meV, Timp=2x 10“^"* s, Tso=10~^^ s, (b) rFM/Tsc=l, ho=30 meV, 
Timp=2xl0“^'‘ s, r.5o=10“^^ s, (c) rFM/Tsc=l, /io=100 meV, rii„p=2x 10“^^ s, Tso=10"^^ s. The 
ferromagnetic and superconducting transition temperatures are denoted by Tfm and Tsc, respec- 
tively. It is interesting that the period increases with T in Fig. 6 (b), while that decreases with T 
in Fig. 6 (c). 
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Fig. 7. The temperature dependence of Ts(T) and he^{T) for Tfm/Tsc = 1 are plotted by the broken 
and solid lines, respectively. 
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Fig. 9. Normalized Josephson critical current as a function of T/Tsc in the case of (a) d=6.5 nm and 
(b) d=7.5 nm, for Tfm/Tsc =17 ^o= 30 meV, Timp=2x s, Tso = 10“^^ s. 
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Fig. 10. The ratios of the period and the decay length, are plotted; (a) by the solid line for 

Tfm/Tsc=10, ho=50 meV, Timp=2x 10“^^ s, Tso=10“^^ s, (b) by the broken line for Tfm/Tsc=1) 
ho=30 meV, Tii„p=2x 10“^^ s, Tso=10~^^ s, and (c) by the dotted broken line for Tfm/Tsc=1) 
ho=100 meV, Timp=2xl0“^^ s, Tso=10“^^ s. The other parameters are set to oJn=o and wf=2.5x 10® 
m/s. The ratio is dominated by T in the case of (a) or around T ~ 0 K in the cases (b) and (c), 
since ttT is larger than 1/ts(T). On the other hand, below Tsc in the cases of (b) and (c), the 
ratio exponentially increases due to T-dependence of 1/ts(T), which is much larger than T. 
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